Normal State Magnetic Properties of Ni and Zn Substituted in YBa2Cu306+x: 

Hole-Doping Dependence 
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We present SQUID susceptibility data on Zn and Ni substituted YBa2Cu306+a;. Cross-checks with 
NMR yield an unprecedented accuracy in the estimate of the magnetic susceptibility associated with 
the substituants, from the underdoped to the lightly overdoped case. This allows us to determine the 
Weiss temperature 6 for YBCO: its value is very small for all hole dopings n^. Since in conventional 
metals, the Kondo temperature, Tk < 0, magnetic screening effects would not be expected for 
T ^ 9; in contrast, increasing Uh produces a reduction of the small moment induced by Zn^"*" and 
a nearly constant effective moment for Ni^"*" corresponding to a spin 1/2 rather than to a spin 1. 



There is widespread interest in the study of substi- 
tutions in correlated magnetic systems. It encompasses 
high Tc cuprates, with the well known depression of Tc, 
as well as cuprate-related materials such as spin ladders 
and even Cu chain materials. All of these systems share 
the features that the local magnetic environment around 
the defect differs from that of the bulk material and cre- 
ates a long range magnetic perturbation. For instance, 
a staggered static magnetization is induced by defects 
in chain compounds H; similarly in cuprates local mo- 
ments appear around Zn substituted in the Cu02 planes, 
as observed by NMR Yet, the underlying theory 
is especially appealing for high Tc cuprates since holes 
play a dominant role and require a different microscopic 
approach. In turn, in this context, the study of sub- 
stitutions in the Cu02 planes helps to achieve a more 
complete understanding of the normal state. 

The differences between various substituants for Cu in 
cuprates, e.g. Zn'^+ (3^1°, 5 = 0) and Ni2+ {3d^, S = 1), 
are still a matter of debate and probably lie not only in 
the magnetic properties of these 2 ions but also in the 
electronic structure of the defect centres, which are re- 
sponsible for the scattering of holes Q . Theories clearly 
suggest that, around the impurities, the underlying an- 
tiferromagnetic correlations of the host are revealed [^. 
Hole-doping dependence is thus expected and is already 
observed in Tc depression ||^ and resistivity 0|. On the 
magnetic side, the focus of this Letter, accurate determi- 
nations of the susceptibility associated with such defects 
have been, up to now, prevented by the existence of spu- 
rious magnetic phases |^,^. Thanks to special care in 
sample preparation, we present here, for the first time, a 
full study of the variation with hole-doping of the effects 
of Ni and Zn substitutions on the magnetism of the cor- 
related Cu02 planes. Our novel experimental determina- 
tion of the Weiss-temperature enables us to demonstrate 
that, for Zn, the weakness of the moment is not related 
to the classical Kondo effect. We present a possible ex- 
planation based on strong correlations. 

Of all the cuprates, YBa2Cu306+2: is the best suited 



for such a study, as the effect of doping on local mag- 
netism can be followed by merely changing the oxy- 
gen content and does not require new synthesis. Our 
YBa2(Cui_yMy)306+a: (M = Zn, Ni) samples, hereafter 
noted Oe+x-M.y%, were prepared at 920°C, using Zn and 
Ni oxides to introduce impurities, with y < 4%. Above 
this limit, substantial amounts (^1%) of spurious phases 
were created. Deoxidizations have been performed under 
thermogravimetric control between 350° C and 430° C in 
a 0.1 mbar vacuum or with a 10 mbar partial pressure of 
oxygen for 0.85 < x < 1.0 (accuracy Ax — 0.01). Magne- 
tization data were taken using a SQUID magnetometer 
in 0.2-5 Tesla fields and Tc's were determined from ac 
susceptibility or SQUID measurements in low field. 

The susceptibility x can be divided into 3 contribu- 
tions, x = Xm + Xybco + Xspur- Thc first one, of primary 
interest here, Xm, is associated with local moments in- 
duced by substitution. It will be shown to follow accu- 
rately a Curie (or Curie- Weiss) law. The second one is 
similar to that of non-substituted YBCO and can be fur- 
ther split into 2 parts, Xybco — XCu{2) + Xch, associated 
respectively with planar Cu far from the substitution site, 
and with CuO chains. NMR, as will be shown hereafter, 
allows one to contrast between XCu{2) and Xch- Finally, 
Curie- like spurious terms due to extrinsic phases, Xspur 
are usually not easily separable from Xm- 

As a prerequisite to the determination of Xm, stud- 
ies of Xybco were first performed on pure powdered or 
sintered samples prepared in the same manner. We suc- 
ceeded in preparing nearly parasitic-phase-free samples 
iXspur ~ 0) as evidenced by the absence of any detectable 
Curie behavior for a; = 1.0 (fig.l). Our accuracy allows a 
refined analysis of Xybcot beyond the (unjustified) usual 
assumption that x is constant. The intrinsic spin sus- 
ceptibility of the planes, measured by ®^Y NMR on our 
sample, decreases by about 7% (~ 2 x 10~® emu/g), from 
150 to 300 K, as commonly observed by NMR |l|l, [||. 
Thus thc observed ~1% quasi-linear increase of Xybcot 
over this T range, fitted as Sxybcoy ~ '^^^'^l^^ from a 
near cancellation of the decreasing planar susceptibility 
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with an increasing chain spin susceptibility (xch)i consis- 
tent with Cu(l) NMR shift measurements The sus- 
ceptibihty of our set of pure samples, all prepared from 
the same master batch, is also reported in the inset of 
fig.l, for 0.55 < X < 1.0. For the most depleted samples 
(x ^ 0.65), Xcu(2) is known from ^''O NMR to decrease 
monotonously at low T Therefore, the flattening 

displayed in the inset of fig.l for x = 0.66,0.55 and T 
< 100 K is not associated with Xspur, as the low-T heat 
treatment performed in order to deoxidize our sample 
cannot create any new parasitic phases. It is the signa- 
ture of a minor paramagnetic contribution due to Xch, 
very likely short chain segments. Finally, for x — 1.0, 
^^Y NMR shift measurements at 100 K on our samples 
show that Zn substitution does not affect the level of hole 
doping in the planes, to within 0.02 in x |16|. We have 
therefore chosen to refer the x value to a; — 1.0 taken 
for the maximum obtained oxygen content, whatever the 
substitution fraction, y. 

The smallness of the value of the induced moment in 
the Zn case makes the analysis of the data far more dif- 
ficult than for Ni. We first study a y = 4% Zn substi- 
tution for X = 1.0 or 0.66. Fig.l shows our data for the 
O7.0 :Znl% and 4% lightly overdoped samples in com- 
parison with the corresponding pure sample. The con- 
tribution of Xm can be extracted by a fit of the T- de- 
pendent part of the raw data to a AT + C/T law where 
the AT term accounts for the contribution of the non- 
magnetic Cu sites (chain + planes) similar to <5Xybco7 ' 
as detailed in the previous paragraph |]l7| . There was 
no significant improvement in the quality of the fit if we 
took a Curie- Weiss C/(T + 0) behavior (-5 < 6* < 10 
K). Fig. 2 shows our results for 06.66:Znl% and 4% un- 
derdoped samples. We focus first on the latter for which 
data could be extended down to 5K, since it is not su- 
perconducting. This makes the Curie increase of x at low 
T more evident. However, the contribution of Xybco (-^) 
to the total T- variation of Xm, is already quite sizeable 
above 60 K (much more than for x = 1.0). Various esti- 
mates of Xybco were therefore performed, 

(i) using the susceptibility of the pure sample, 
XYBcogeje' Pi'^sented in the inset of fig.l (dotted line), 
which limits the fit to the T-range T >75 K 

(ii) using the data from planar ^^O or ^^Y NMR mea- 
surements rescaled by the hyperfine constant and extrap- 
olated linearly to zero spin shift at T = 0. 

The main advantage of (i) is to include contributions 
from the chains, (xc/i), at low T |l^ (see fig.l) whereas 
in (ii) Xch is completely neglected. The local moment 
contribution, Xm X ^ (1 ^ 5c)Xybco is then deduced, 
yielding the correct Curie- Weiss fit over the whole T- 
range. The factor 1 — 5c (c = 1.5y is the concentration 
per plane) reflects the fact that the substituted Cu site 
and its four neighbors which carry the local moment are 
not expected to have the same susceptibility Xybcos bJ-?^ 
the bulk Cu of the CuOa planes (as shown by Y NMr'[§). 
The high-T fit (procedure (i)) yields C = 2.3(2) x 10^^ 
emu.K/g and 6* = 5(3) K whereas, including the low-T 



data gives (ii) C = 2.5(2) x IQ-'^ emu.K/g and = 4(1) 
K. C is modified by less than 15% if we extract Xm by 
removing the contribution of the substituted site only. 



that is (1 — c)x-y 



The dominance of the Zn-induced 



contribution at low T (T < 30 K) makes the value of 6 
hardly dependent on the T = extrapolation performed 
for Xybco i^ approach (ii) . It is remarkable that all meth- 
ods yield very similar results within a 15% error bar for 
a given sample. The value of 6 was further checked by 
studying the variation of the magnetization M with the 
applied field H for 5 < T < 12 K (fig. 3); indeed, one ex- 
pects M to saturate for high fields and low T, according 
to the Brillouin function M = Bj{^iBH/kB{T + 6)). We 
find 9 — 4.5(5) K, in good agreement with the previous 
values. Finally, it is interesting to note that 9 is quite 
small and comparable to the transition temperature of 
the disordered magnetic state (~2 K) observed in fiSK 
for the same 06.66:Zn4% sample Q. 

Determinations of C were also made along the same 
lines for lower Zn substitution rates but the occurrence 
of superconductivity prevented us from estimating 9 as 
accurately as for y = 4%. In any case the data plotted in 
fig. 4 clearly demonstrates that 9 never exceeds a few K. 

Our complete experimental data analysis for Zn and 
Ni, partially reported on fig.l, 2 and 4, is summarized in 
the inset of fig. 4. We find that C increases linearly with 
y. The values of the local moment are deduced from a 
linear fit of C(y) ~ Ul^lff and results are given in Table 
I, along with the depression of Tc- One of our important 
experimental conclusions is that the value of the Ni mo- 
ment stays quasi-constant with hole doping and near a 
spin 1/2 value (/^e// = 1 -73/15), far from the expected 
5 = 1 for a Ni^+ ion (/ie// — 2.82/iB). In contrast, the 
moment induced by Zn is quite small. C decreases by 
a factor 5.8 from the under- to the lightly over-doped 
side of the YBCO phase diagram. One might wonder 
whether the reduction of the moment for Zn is due to 
a classical Kondo effect. In ordinary metals, such a re- 
duction only occurs below 9, which is of the order of 
the Kondo temperature Tk- On the contrary, here, for 
T ^ 9, the moment is reduced whereas Xm deviates only 
slightly from a Curie law (moreover 9 is rather associated 
with long-range interactions between moments). 

Finally, we monitored the decrease of the Curie con- 
stant, C, with hole doping for the 4% Zn concentration. 
In the upper panel of fig. 5, we evidence a strong corre- 
lation between C and the residual resistivity po taken 
from 10]. This experimental finding suggests that local 
moments and scattering are closely related for cuprates, 
which, again, is at variance with the case of ordinary 
metals where non-magnetic impurities induce scattering 
without creation of any related magnetism. Accordingly, 
models should account not only for the scattering cross- 
section but also for the value of the moment induced by 
the substituting element. We next present a tentative 
theoretical interpretation, based on strong correlations, 
which relates these two quantities. 

First, we consider the Zn case where the impurity can 
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be treated as a spin vacancy. Simulations for a weakly 
doped t-J model for cuprates yield a local moment for 
Zn 1^, which agrees qualitatively with our findings but, 
for our experiments, the hole content is far larger. The 
extension to higher doping can be treated in the frame- 
work of gauge theories which advocate spinon-holon sep- 
aration. Using this approach, Nagaosa and Lee po|[_w ere 
able to account for resistivity data on 06+a::Zii|7[], in 
the two extreme cases = (overdoping) and 5=1/2 
(underdoping) . Following Gabay jl^, a description cov- 
ering the whole doping range can be obtained. Because 
Zn is in a Sc?^" configuration, a singlet is always present 
on the substituted site. According to the t-J picture, this 
implies that both holons and spinous are scattered off 
this site. Their density is thus depleted within a typi- 
cal radius, i?s, around Zn, or in mean field theory within 
Atf, the Thomas- Fermi length; thus the spin 1/2 carried 
by those Cu surrounding the Zn impurity are not fully 
screened and a local moment S survives in the vicinity 
of Zn. For a hole fraction S, the spinon number per unit 
cell is 1 — (5, far from the Zn site. Within a disk of radius 
a (the lattice constant) centered on Zn, it is reduced to 
(1 — S)ps{a) with, for a featureless Fermi surface. 
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with NMR experiments The magnitude of the mo- 
ment is then expected to be independent of doping. In- 
terplane coupling involving the 1^322 orbital may also 
help explain the experimental observation of a small vari- 
ation of S with hole doping. 

To conclude, we have found a strong contrast between 
the case of Ni and Zn substitutions in YBCO. Not only 
the depression of Tc and the value of S are drastically dif- 
ferent but also the interplay of substitution with charge 
doping is now obvious in the case of Zn. This sets strong 
constraints on theoretical models by requiring them to 
link local moments and scattering by impurities, beyond 
the classical Kondo picture. We have addressed this point 
using a model based on correlations but the problem is 
still open as, e.g., the variation of the magnetic corre- 
lation length with doping could also play a role. De- 
viations from a Curie law are found negligible and more 
likely due to a very small spin freezing temperature linked 
with interaction between moments rather than classical 
Kondo behavior. The smallness of the moment observed 
in OriZn raises the fundamental question: does it persist 
into the heavily overdoped phase of the cuprates? 

We thank P. Monod and S. Zagoulaev for discussions 
and preliminary ESR tests of our sample quality. 
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For scattering states in 2D characterized by a wavevec- 
tor fcoj Rs is related to the planar resistivity, po,2D by 
Rs = e^6/2koh x po^2d/c. In the gauge picture, the 
main contribution to po.2D comes from holons, that is, 
neglecting band effects, fco = -\/27r(5/(?. Connecting 
Rs to Xtf gives the variation of the effective moment 
Peff = gi^J'B\/S{S + 1) with respect to doping. Stan- 
dard scattering theory, when applied in 2D, gives that 
Rs/Xtf ~ S-^^^ pa so that 
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From the data of ref. |Q for Po,2d/c (expressed in 
m/D.%Zn) and Eq||, we find that a =0.050(5) fits sat- 
isfactorily our data for fieff{x), in fig. 5. A more refined 
treatment [beyond the scope of this Letter] would include 
taking into account the actual shape of the Fermi surface. 

The case of a magnetic substitution (Ni) is quite dif- 
ferent: no magnetic vacancy is introduced into the Cu02 
plane. A theoretical approach paralleling that for Zn 
impurity scattering also accounts for resistivity measure- 
ments in the case of Ni doping. In this latter case, scat- 
tering is due to the formation of a singlet bound state on 
the Ni site. In addition there remains a 3 d^z^-r'^ S = 1/2 
spin on Ni, hardly coupled to the planes, in agreement 
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07:Ni 06.66:Zn 


06.66:Ni 


(K/%) 


11.0(1.5) 


3.4(7) 21(2) 


7.5(1) 




0.40 (5) 


1.2(1) 1.0(1) 


1.6(1) 



TABLE I. Tc depression and local moment for Zn and Ni 
substitutions in lightly overdoped and underdoped YBCO 



FIG. 1. Susceptibility of 07:Znl% and 4% samples com- 
pared to pure 07(open squares). Notice, for the latter, the 
increase of x with T; whereas, the susceptibility of the Cu02 
planes (Y NMR, shifted up by lO"'^ emu/g) decreases (dashed 
line). Inset Xybco6+, Jox x = 0.56, 0.65, 0.76, 0.85, 1.0. 



FIG. 2. Lower panel: High-T plot of x, Xm for 06.66:Znl% 
and 4%. The closed symbols stand for raw data (x), the open 
ones for Xm, obtained after subtraction of Xybco using (i) 
(see text). The lines are Curie- Weiss fits described in the 
text. Upper panel: Low-T plot of Xm(YBC06.66:Zn4%) vs. 
1/{T + 6) {6 = 4.0 K). The line corresponds to the same value 
of C as determined in the lower panel. 

FIG. 3. Plot of the magnetization versus H/ (T -I- 9) for 
various values of 9. The upper and bottom panels show the 
deviations from the expected universality for a Brillouin func- 
tion, for unappropriate values of ^ (^ = 2 and 7 K). The uni- 
versality is only obtained for a narrow range of values of 9, 
9 = 4.5(5)K (mid-panel). 

FIG. 4. Xm versus 1/T for 06+i:Zny%,Nij/%. The lines are 
fits described in the text. Inset: Curie constant (in emu.K/g) 
versus (%) impurity content for x = 7.0 and x — 6.66. 

FIG. 5. Top panel: Evolution of C ior y = 4% (closed 
squares), and po,2D (open squares) with oxygen content. 
Bottom panel: Evolution of /ie/ / (closed squares) and Tc (tri- 
angles), for y — 4%. The open squares are estimates for /leff 
using eq. | for S. 
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